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Abstract. Using time- resolved two-dimensional aperture 
photometry we have put upper limits on the pulsed emis- 
sion from two proposed optical counterparts for PSR 
B1951-I-32. Our pulsed upper limits of m„ puised > 23.3, 
'Tife.puised > 24.4, for the first candidate and m„_puiscd > 
23.6, TO;,, pulsed > 24.3 for the second, make it unlikely that 
either of these is, in fact, the pulsar. We discuss three fur- 
ther candidates, but also reject these on the basis of tim- 
ing results. A search of a 5'.' 5 x 5" 5 area centred close to 
these stars failed to find any significant pulsations at the 
reported pulsar period. 

Key words: pulsars: individual: PSR B1951+32 



1. Introduction 

PSR B1951-I-32 is a fast pulsar in the peculiar combina- 
tion supernova remnant (SNR) CTB 80. The characteris- 
tic spin-down age of the pulsar {P/2P) and the dynam- 
ical age of the SNR (Koo et al. 1990) are both - 10^ 
yrs. Its low surface magnetic field strength, of 5 x 10^ 
T, is unusual amongst the selection of young and middle- 
aged pulsars detected so far but it has been suggested 
that this is more representative of young pulsars in gen- 
eral ( Lyne et al. 1983| ). It is clearly important to test pul- 
sar models with as wide a variety of pulsar parameters as 
possible. 

Steady emissi on from the pulsar has so far bee n de- 
tected at radio ( ^trom 1987 ; Kulkarni et al. 1988 ) and 
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X-ray ( [Becker et al. 1982[ [Wang fc Seward 1984D wave- 
lengths. Two possible optical counterparts have been 
proposed by Blair & Schild (1985) and Fesen & Gull 
( 1985). Since the ini tial discovery of the 39.5-ms pulsar 
( Kulkarni et al. 1988 ), evid ence for pulsed emis sion has 
also been found in X-rays ( ^afi-Harb et al. 1995 ) and 7- 
rays (Ramanamurthy et al. 1995) with upper limits in the 
infrared (Clifton et al. 1988). To-date no upper limits on 
pulsed emission in the optical have been published. 

Optical observations are essential in determining the 
relative contributions of thermal and magnetospheric pro- 
cesses to pulsar emission. The continuing improvements in 
optical-detector sensitivities means that pulsed emission 
from 6 isolated neutron stars has already been detected. In 
this paper we present time-resolved optical observations of 
the central region of the supernova remnant CTB 80 and, 
in particular, of the two proposed pulsar candidates. 



2. Observations 

The observations were made during 1996 June and 1997 
June using NUI, Galway's TRIFFID camera. The TRIF- 
FID system consists of a multianode microchannel array 
(MAMA) 2-dimensional photon-counting detector with a 
B extended S-20 photocathode ( flimothy fc Bybee 1985| ) 
and a fast data-collection system ( pledfern et al. 1993 ). 
The position and time-of-arrival of each photon are 
recorded to a precision of 25/im and 1/is, respectively. 
TRIFFID's 1024 X 256-pixel array was mounted at the 
prime focus of the SAO 6-m Telescope (BTA), resulting 
in an equivalent spatial resolution of 0'.'23 pixel" ^. Timing 
was achieved using a GPS receiver and an ovened lO-MHz 
crystal. Absolute timing was unavailable for the 1996 ob- 
servations however, and times were recorded relative to 
the last reset of the data-collection system. The data- 
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collection system was reset at the start of the observing 
run and again after the 1996 June 23 observations. 

The observations were made with standard B and V 
filters and are summarised in Table 1. In all, a total of 
19,456 seconds of B-band data and 8,173 seconds of V- 
band data were collected. Photometric standard stars were 
observed each night for calibration purposes. 



Table 2. PSR B1951+32 Ephemeris (Foster et al. 1994) 



Parameter 


Value 


V 


25.29739673Hz 


V 


-3.73571 X 10"^^ Hz 


V 


-3.37 X 10"^^ Hz 


Epoch 


2447005.6880 



Table 1. Summary of Observations 



Date 


UTC 


Filter 


Duration 


Seeing 








(s) 


(") 


1996 June 22 


20:31:00 


B 


9013 


1.40 


1996 June 23 


21:23:00 


V 


3756 


1.16 


1996 June 23 


22:28:30 


B 


1333 


1.28 


1996 June 23 


22:52:40 


B 


3110 


1.45 


1996 June 24 


22:33:41 


V 


4417 


1.40 


1996 June 25 


23:00:21 


B 


3300 


1.24 


1997 June 04 


22:59:17 


B 


2700 


1.50 



3. Data Reduction 



The data were first binned into 1-ms frames and di- 
vided by a deep flatfield image taken during the obser- 
vation period. A post-exposure shift-and-add sharpening 
technique was applied to produce an integrated image 
( ^hearer et al. 1996 ). Because of the large telescope aper- 
ture no significant improvement in the image above the 
seeing limit resulted, but any artefacts due to effects such 
as telescope wobble were removed. Data from individual 
nights were summed to produced the integrated V and B 
images shown in Figs. 1 and 2. 

Photometry and astrometry were carried out using 
the IRAF daophot and GASP packages, respectively. The 
background level was taken to be the mean of the signal 
in an annulus of radius 2'.'0 and width O'.'S centred on the 
object position. Stellar co-ordinates were calculated using 
the astrometry of Blair & Schild (1985). 

Several pulsar candidate positions were chosen as de- 
scribed in Sec. 4. Photon times were extracted from 
a window with a diameter equal to the seeing width 
centred on the pulsar candidate. The time series was 
translated to the solar-system barycentre using the JPL 
DE200 ephemeris and t hen folded in phase using the PSR 
B195H-32 ephemeris of [Foster et al. 1994| (Table 2). The 
resulting light curves were analysed using the Z"^ statistic 
( Buccheri &: de Jager 198E ). The Fourier power-spectrum 
of each candidate was also calculated in the vicinity of the 
first four harmonics. 



4. Results 



Figs. 1 and 2 show the integrated V- and B-band images. 



The proposed opti cal counterparts ( Blair fc Schild 1985 ; 
Fesen & Gull 1985) labelled '1' and '2', can clearly be seen 
in both images. The core region surrounding the two can- 
didates is shown in Fig. 3, which is a sum of all the observa- 
tions. The positions of the radio pulsar (Foster et al. 1994) 
and X-ray point source (Wang & Seward 1984) are also 
marked. Like Fesen & Gull we find no evidence for Blair 
and Schild's candidate '3' (down to a limiting magni- 
tude of m-u ~ 24), consistent with its identification as a 
Ha -I- [Nil] emission knot. The photometry for these can- 
didates, and those discussed later on, is given in Table 3. 
The magnitudes were calibrated with respect to Blair and 
Schild's star 'A'. The quoted errors were calculated from a 
combination of counting statistics and flat-fielding errors. 
Our B1950 source positions of the two pulsar candidates 
are a = 19''51"02!58, 5 = -|-32°44'50'.'9 (candidate 1) and 
a = 19''51™02!49, (5 = -|-32°44'52'.'0 (candidate 2). 

Estimates of the distance to the pulsar vary between 
about 1.4 kpc (from dispersion measure estimates, Kulka- 
rni et al. 1988) to 3 kpc (studies of the infrared shell, Fesen 
et al. 1988). Following Safi-Harb et al. (1995) we adopt 2.5 
kpc as a compromise between these extremes. Reddening 
estimates also vary and lie within the range E(B-V) = 
0.8-1.4 (Angerhofer et al. 1980; Blair et al. 1984). We use 
a value E(B-V)= 1.0 which also gives a distance to the 



pulsar of about 2.5 kpc (Blair & Schild 1985). The mea- 
sured values are equivalent to a star of absolute magnitude 
Mt, — 5.30 and intrinsic colour (B-V)o — 0.77 (candidate 
1) and M, = 6.29, (B-V)o = 0.86 (candidate 2). 

Photon times were extracted from apertures centred 
on the two pulsar candidates. No pulsations were seen, 
to the 1% significance level, in either the B- or V-band 
data sets. An upper limit for the pulsed fraction was cal- 
culated for the B and V emission assuming a duty cycle in 
the optical of 50%. These upper limits are conservative as 



pulse profiles in the X-ray (Safi-Harb et al. 1995) suggest 



a much lower duty cycle. The 3- a upper limits were found 
to be TO^^puisod > 23.3 ± 0.41, TOfc^puiscd > 24.2 ± 0.35, for 
candidate 1, corresponding to pulsed fractions < 6% in 
V and < 11% in B. For candidate 2 we found m^^puised > 
23.6±0.41, mfj^puisod > 24. 2±0. 35, corresponding to pulsed 
fractions < 12% in V and < 32% in B. These limits were 
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Table 3. Photometry for bright sources and pulsar candidates in the CTB 80 field. (The values for star A are taken from Blair 
& Schild (1985)). Vahios for M,. fc (B-V),, ar(> c-laculatod asuiriiiig' a distance to the ol)j(X't of 2.5 kpc. 





R.A.(1950) 


Dec. (1950) 












Source 


h m s 


/ // 


V 


B 


(B-V) 




(B-V)o 


A 


19 51 02.10 


32 45 15.0 


14.95 ± 0.20 


16.86 ± 0.20 


1.91 






1 


19 51 02.57 


32 44 50.42 


20.29 ± 0.30 


22.06 ± 0.20 


1.77 ±0.30 


5.30 


0.77 


2 


19 51 02.47 


32 44 51.83 


21.28 ± 0.30 


23.14 ±0.30 


1.86 ± 0.30 


6.29 


0.86 


4 


19 51 02.55 


32 44 49.25 


22.08 ± 0.20 


22.87 ±0.20 


0.79 ± 0.30 


7.09 


-0.21 


5 


19 51 02.65 


32 44 48.92 


22.55 ±0.30 


23.12 ±0.20 


0.57 ±0.30 


7.56 


-0.43 


6 


19 51 02.48 


32 44 47.63 


23.05 ± 0.30 


23.01 ± 0.20 


-0.04 ± 0.30 


8.06 


-1.04 




Fig. 1. V-band image of the core region of CTB 80. The image 
is a mosaic of two separate observations lasting 8,173 seconds 
in total. The star identifications refer to those listed in Blair 
& Schild (1985). (1 pixel = 0'.'23) 



calculated from the longest series of data taken with a 
given filter between resets of the data-collection system 

(sec Section 2), a single data file (4,417 seconds) in the 
case of V and the sum of three files (13,456 seconds) in 
the case of B. 

In Fig. 3 candidate 1 appears to be extended in the 
direction of the radio pulsar position, a feature that can 
be seen in both V and B data sets separately. 

An accurate PSF for the image was determined us- 
ing the external IRAF package DAOPHOTX, and used to 
remove the candidate 1 and 2 stars from the image. We 
found that the extension was consistent with two point 
sources. Fig. 4 shows the extension along with a another 



Fig. 2. B-band image of the core region of CTB 80. The image 
was made from the B-band observations listed in Table 1, and 

is equivalent to a total exposure of 16,456 seconds. The star 
identifications refer to those listed in Blair & Schild (1985). (1 
pixel = 0'.'23) 



point-like feature nearby (we have labelled these candi- 
dates '4','5' and '6' (Fig. 4). 

We have tried to assign a magnitude and colour to each 

of these features, although their low signal-to-noisc ratio 
and proximity to a brighter source makes this difficult to 
do with any accuracy. The photometry for these sources 
is included in Table 3. Absolute magnitudes and intrinsic 
colours are calculated assuming a distance of 2.5 kpc to 
the pulsar. 

Times were extracted from all these candidates and 
analysed using the Z'^ statistic and their Fourier power 
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Fig. 3. V+B image of the central region of CTB 80. Blair and 
Schild's candidates 1 and 2 can be seen toward the centre of the 
image. The radio (white cross) and X-ray (black cross) source 
positions are marked. (1 pixel = 0'.'23) 



spectrum. Again, no pulsations were seen to the 1% sig- 
nificance level. 

Fesen & Gull's (1985) Ha-f [Nil] map of the core shows 
that some nebulosity is coincident with the candidate po- 
sitions, although little was detected by Blair & Schild 
(1985). Spectra of CTB 80 filaments ( |Blair et al. 1984] ) 
show strong [01], Ha, [Nil] and [SII] emission lines, but 
these are removed by our filters. Emission lines in the B 
and V bands are much weaker, although our images do 
show some evidence of nebulosity. We cannot rule out the 
possibility that these fainter point-like sources are emis- 
sion knots in the nebula. 

Finally, we extracted photon times from apertures of 
radius O'.'Sl centred on a two-dimensional grid of points 
separated by 0'.'46 on the 1997 data. The grid covered an 
area 5'.'5 x 5'.'5 centred close to the first candidate posi- 
tion. The photon times were reduced to the solar-system 
barycentre as before, and the Zf statistic calculated for 
each aperture using the ephemeris in Table 2. No evidence 
for the pulsar was found. 
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Fig. 4. Contour plot of the region surrounding candidate 1 in 
the V-band image of Fig.l. The position of the radio pulsar 
(Foster et al. 1994) is marked with a cross. (1 pixel = 0'.'23) 



the neutron star surface itself, as we estimate that a pulsar 
with T « 10^ K at the distance of CTB 80 would have a 
magnitude in the range mi,=30-31. We may, however, be 
able to see emission from any hot circumstellar material. 

From phenomenolgical models of pulsar optical emis- 
sion (Pacini and Salvati 1983, 1987) wc can derive esti- 
mated magnitudes for the emission. For PSR 1951+32 
this leads to values in the range 24-26. However recent 
observations of other middle-aged pulsars, 0656-1-14 and 
Geminga (Caraveo et al. 1994, Shearer et al. 199S), have 
indicated a magnitude considerbly in excess of that pre- 
dicted by the Pacini and Salvati model. It is clearly very 
important to establish a positive optical identification for 
PSR B1951-I-32 in order to compare emission from pul- 
sars of similar ages but with a range of other parameters 
such as P and magnetic field. These observations will be 
crucial if we are to distinguish between the various mod- 
els for high-energy emission. In particular, the pulse shape 
and fraction can be used to differentiate the polar cap and 
outer gap models. 

We suggest observations with the VLT, Keck or the 
HST in order to definitively identify the optical counter- 
part of PSR B1951+32. 



5. Conclusions 

On the basis of luminosity, colour and timing we conclude 
that neither of the proposed optical counterparts are the 
pulsar PSR B1951-h32. 

We note that the slight blue extension to candidate 
1 shows the expected characteristics of a hot object. We 
would not expect to be sensitive to thermal radiation from 
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